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The purpose of this paper is to examine the possible existence or construction of traversable worm- 
holes supported by generalized Chaplygin gas (GCG) by starting with a general line element and 
the Einstein tensor, together with the equation of state, thereby continuing an earlier study by the 
author of wormholes supported by phantom energy. Numerical techniques are used to demonstrate 
the existence of wormhole spacetimes that (1) meet the flare-out conditions at the throat, (2) are 
traversable by humanoid travelers, thanks to low tidal forces and short proper distances near the 
throat, and (3) are asymptotically flat. There appears to be an abundance of solutions that avoid 
an event horizon, suggesting the possibility of naturally occurring wormholes. 
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I. INTRODUCTION 

Interest in traversable wormholes, initiated by Morris and 
Thorne [1] , has been greatly renewed in part because of 
the discovery that our Universe is undergoing an accel- 
erated expansion [3, 's']. In other words, a > in the 
Friedmann equation a/ a — —^{p + 3p), using units in 
which G ~ c=l. The cause of this acceleration is taken 
to be a negative pressure dark energy with equation of 
state p = —Kp, K > ^, and p > 0, where K is a. con- 
stant, p the spatially homogeneous pressure, and p the 
energy density. Of particular interest is the case K > 1, 
referred to as phantom energy, since it leads to a vio- 
lation of the null energy condition, an essential require- 
ment for maintaining a wormhole [ij. Matter violating 
the null energy condition is usually called exotic. Since 
the notion of dark or phantom energy ordinarily applies 
only to a homogeneous distribution of matter, phantom 
energy is not automatically a candidate for exotic mat- 
ter. Fortunately, the extension to spherically symmetric 
homogeneous spacetimes has been carried out [4]. 

An alternative model is based on Chaplygin gas, whose 
equation of state is given hy p = — ^ . Another possibility 
is generalized Chaplygin gas (GCG), whose equation of 
state isp=— ^,0<a<l Cosmologists became 

interested in this form of matter when it turned out to 
be a candidate for unifying dark matter and dark energy. 
To see this, consider the energy conservation equation 
jj — —3d{p+p)/a in a flat FRW spacetime and substitute 
the equation of state p = The result is 

where _B is a constant of integration. It is now seen 
that p ^ at early times, that is, p behaves like 
matter, while in later times like a cosmological constant 



(p — constant). Previous models required two distinct 
fields, one to describe dark matter and the other dark en- 
ergy, but one can argue, as in Ref. 0, that these ought 
to be different manifestations of the same entity. One 
possible motivation for this model came from the field 
theory points of view l8jj . Another interesting possibility 
is discussed in Ref. Q: starting with the Nambu-Goto 
action of string theory, the Chaplygin gas appears after 
considering c?-branes in a (d -I- 2)-dimensional spacetime. 
Another attraction is that it admits a supersymmetric 
extension. 

Obscrvationally, the GCG model has not been with- 
out its problems. As noted in Ref. [l3|, while the 
model has successfully withstood various phenomenolog- 
ical tests over several years, there is some concern that 
it produces unphysical oscillations in the matter power 
spectrum. It is shown in Ref. [13], however, that these 
problems can be circumvented. Furthermore, very recent 
studies [ill . IT^ have concluded that the earlier criticisms 
were based on the oversimplifying assumption of an adi- 
abatic cosmic medium. 

On the other hand, according to another recent study 
(Ref. [l3|), the special case a = 1, corresponding to 
the original Chaplygin gas, may very well have to be 
excluded. One should therefore concentrate on the GCG 
case by keeping a < 1. 

In a recent paper [l^l the author made a systematic 
study of exact solutions of wormhole spacetimes sup- 
ported by phantom energy by starting with the general 
line element and equation of state. It is shown that there 
are only two ways to insert the redshift function "by 
hand." Doing so leads to the exact solutions in Refs. [l5j 
and 16]. Assigning a specific function to p leads to the 
exact solution in Ref. [Tj]- Additional exact solutions 
that simultaneously avoid an event horizon are extremely 
rare. Nevertheless two new solutions were found. In- 
cluded in the discussion are the junction conditions for 
matching each solution to an exterior Schwarzschild so- 
lution, as well as traversability criteria. 

Ref. 19] makes use of specific redshift and shape func- 
tions in the manner of the original Morris and Thorne pa- 
per. In this paper we continue in the spirit of Ref. [14| . 
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that is, we start with the general hne element and the 
Einstein tensor, together with the equation of state, and 
continue the analysis while avoiding, as much as possible, 
assigning specific functions to the metric coefficients. 



II. THE PROBLEM 

Consider the general line element for describing a worm- 
hole: 



2 + e2A('^)rf 



+r2(d6|2 +sin2 6* #2). (2) 



In this context $(r) is called the redshift function; for 
this function, e^**^''^ must never vanish to avoid an event 
horizon. A(r) is related to the shape function b{r) = 
r(l - e-2AW)^ i.e., e^^^ = 1/[1 - b{r)/r]. The shape 
function determines the spatial shape of the wormhole as 
viewed, for example, in an embedding diagram. By the 
very definition of wormhole, if the throat is at r = tq, 
then b{ro) — tq. As a consequence, A(r) has a vertical 
asymptote at r = rg: limr^ro+ A(r) = +oo. To obtain a 
traversable wormhole, the shape function must obey the 
usual flare-out conditions at the throat 
and b{r) < r; also required is asymptotic flatness, i.e., 
h{r)/r ^ as r ^ oo. 

The next step is to list the components of the Einstein 
tensor in the orthonormal frame Il4|: 

Gu - 'e--(^)A'M + 1 (l 



-2A(r) 



(3) 
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-2A(r) 
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g-2A(r) _ $'(r)A'(r) + [^'{r)f 



+ i$'(r) - -A'(r) ) . (5) 
r r J 

Since the Einstein field equations — SttT^^ im- 

ply that the stress-energy tensor T^^ is proportional to 
the Einstein tensor, the only nonzero components are 
Tit = Pi^)^ = p{r), and Tg^ = T^^ = pt{r), the trans- 
verse pressure. From the Einstein field equations and 
the equation of state p = we obtain G^^ — Snp and 

Gff ~ 8tt{—Ap~°'), yielding the following equation: 



1 

8^ 



-e-2AW$'(r)- J-(l 



,-2A(r) 



A 



{«^[^-^^WA'(r) + ;L(i_e-2AW)]}- 



(6) 



Unlike the phantom energy case (p = —Kp), the equation 
of state has the form of a quotient. The resulting differ- 
ential equation ^ does not have an obvious exact solu- 
tion, so that the analysis depends on numerical/graphical 
techniques. 



Recall next the null energy condition, which requires 
the stress-energy tensor Tafi to obey Taf^fi^p,^ > 
for all null vectors. In our orthonormal frame, for 
(/i*, /i'', 0, 0) = (1,1,0,0), a radial outgoing null vector, 
the condition becomes Tii + Tff = p + p > 0. Wormholes 
must necessarily violate this condition at the throat [l|. 
Since limr^ro+ A(r) = +oo, Eq. ^ now implies that 



A 



= 0. 



Solving for A, we have A — l/(87rrQ)"+^. We will see 
below that we actually need 



A < 



1 



(87rr2)«+i • 
Using the above equation of state, 
A 

0<a<l, 



(7) 



(8) 



in conjunction with p + p < 0, yields the constraint 

p<Ai/("+i\ (9) 

which, in turn, implies that _B < in Eq ([T]). (The 
constraints ([7]) and ^ are also discussed in Ref. [19|.) 

III. THEORETICAL ANALYSIS 

As already noted, Eq. ^ does not lend itself to finding a 
simple exact solution, so that some numerical techniques 
will be needed. That is the topic of Sec. IIVI In this sec- 
tion we determine some general characteristics of GCG 
wormholes. The analysis is based on the assumption that 
a solution A(r) of Eq. ^ exists for a wide range of choices 
for <I>'(r) and such that limr^ro+ A(r) — +oo. We also 
assume that $'(r) is continuous. 

Suppose we rewrite Eq. ^ as follows: 

(87r)i+i/°r(r^)i/°Ai/° 

(1 _ e-2A(r-)) [i _ e-2A(r) „ 2re-2A W $' (r)] 

_ 1 2e-^^('-)A^(r) 
~ 1 - e-2A(r) • 

Now define the dimensionless function 



(10) 



Fir) 



(87r)i+V"r'[(r')2]i/"Ai/«dr' 



ro (1 - e-2A('-')) [1 - e-2A('-') - 2r'e-2A('-') $' (r')] ' 

(11) 

Observe that since e"^^^''^ ^ as r — > ro+, the inte- 
grand is sectionally continuous. This implies that F(r) 
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is defined and continuous on any interval (ro,r]. More- 
over, F(ro) — 0. 

Returning to Eq. (fTO|) . integration yields 

F{r) = In (l - e-2A('-)'j + i^j, + i^c. 



So 



(12) 



Once again, liin,.^ro+ ^k^") = +oo ^-nd F^vq) = 0. So by 
Eq. inn, c = l/ro and 



„2A(r) ^ 



1 - la e^('~) ' 

r 

The line element then takes on the form 



(13) 



1 _ m (,F(r) 



The next step is to show that the flare-out condi- 
tions are satisfied at the throat. Given that 6(r) = 
r (l — e"^'^^'')), it follows that b{r) < r for r > tq. Also, 
Eq. ([131) implies that 



b(r) = r f 1 - 1 + ^e^('-) 



and 



bir) = roe^M. 



By Eq. 1^, 
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1 
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then 



< 



(87r)i+i/"(r2)i+i/"' 



whence b'{ro) < 1. 

In summary, the analysis in this section depends pri- 
marily on the asymptotic behavior of A(r) at the throat. 
It has not been established that the existence of A(r) 
yields a traversable wormhole. In addition, its depen- 
dence on $(r) has to be closely examined, just as it is for 
phantom energy supported wormholes |14l |. These ideas 
are discussed in the next section. 



IV. NUMERICAL ANALYSIS 

In this section we turn to the solution of Eq. (fTO|). rewrit- 
ten as follows: 

A'(r) = 

i(87r)i+i/"r(r2)i/«Ai/"e2^('-) 1 / ^^^^^ ^ 
[l - e-2A('-) - 2re-2A('-)$'(r)] 2r T 

(14) 

Here A is taken to be less than l/(87rrQ)"+-'^ and a to be 
close to unity. (In fact, for convenience we let a = 1, at 
least for now.) In order to generate numerical/graphical 
output, it is necessary to choose some specific value for 
the size of the throat (such as ro = 2) and to pick an ini- 
tial value [such as (2.000001,5)] to simulate the asymp- 
totic behavior. As one might expect, qualitatively, the 
form of A(r) remains roughly the same for a large range 
of $'s. 



A. Regression and traversability 

In this section we turn to the question of traversability 
by humanoid travelers. One concern is the calculation of 
the proper distance from the throat of the wormhole to 
a point outside: 



(15) 



For this calculation a specific A(r) would be desirable. 
Regression equations from the numerical output are easy 
to obtain but tend to produce a poor fit — with one no- 
table exception: A(r) = a + bhi(r — ro). This form yields 
an excellent fit for any $(r) with small |<I''(r)| and an al- 
most perfect fit for $'(r) = 0. {^'{r) can be either posi- 
tive or negative; as |$'(7')| increases, A(r) = a+b\n{r—ro) 
yields an ever poorer fit.) Some of the regression equa- 
tions using To = 2 with the corresponding <&' are given 
next: 



$'(r) = : 


A(r) 


= 1.20 


-0.231n(r 


-2) 




A(r) 


= 1.25 


-0.221n(r 


-2) 


$'(r) = lOe"'' : 


A(r) 


= 1.32 


- 0.20 ln(r 


-2) 


= 4 


A(r) 


= 1.45 


-0.181n(r 


-2) 




A(r) 


= 1.57 


-0.161n(r 


-2) 



Observe that each A(r) has a vertical asymptote at r — 2 
and that A'(r) < for r > 2. The numerical solutions 
were obtained only in the vicinity of the throat. The 
reason for this is discussed in Subsec. C. 

Remark 1. We may assume that the equations are 
fairly typical for GCG wormholes since there is nothing 
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special about the choice of rg or the initial conditions. 
Also implied is that the wormholes are macroscopic. 

From the standpoint of traversability, a relatively small 
|$'(r)| results in relatively small tidal forces [ij. In fact, 
for $'(r) = 0, the radial tidal force is zero. 

The proper distance from the throat r = tq to a point 
outside, say r = kro for some fc > 0, is also going to be 
relatively small thanks to the resulting In-form. As an 
example, for A(r) = 1.20 — 0.23 ln(r — rp), we have from 
Eq. lEl) 



fcro 



g..20-0.231n(r-ro)^^ 



kro 



0.77 ' 

(16) 



corresponding to the coordinate distance [tq , kro] ■ 

We conclude that a wormhole with small in the 

vicinity of the throat is likely to be traversable in the 
sense of having low tidal forces and short proper dis- 
tances resulting from the In-form. A wormhole is there- 
fore traversable in this sense if, and only if, the numerical 
output for A(r) closely fits the form A(r) = a-|-61n(r— rg). 

Remark 2. Although of primary importance, |^''(?')| is 
not the only concern involving the redshift function. Re- 
turning to line element ([2]), even if limr_^+oo ^(r) = 0, a 
large near the throat is also undesirable: depend- 

ing on the sign of clocks fixed at r = rp will either 

run much faster or much slower than clocks outside the 
wormhole. So our definition of traversability does not 
include a study of proper traversal times. 



Suppose we consider the expansion 



ln(r 



^o) = lnr-^--:|-i:|- 
2 3 r-^ 



(17) 



obtained by expanding (c?/(ir)ln(r — rg) = l/(r — rg) in 
a geometric series and then integrating. Returning to 
Eq. (HID, 



i{r) = 



kro 



A.20g-0.23[lnr-Er^(^)1d. 



(18) 



which is equal to the value of £{r) in Eq. (|16p . To measure 
the deviation from this form, let us call A(r) "log-like" if 
its expansion is 



r 2P r2 3^ 



(The condition p < 1 makes the function "less favorable" 
than the In- function.) The series itself converges for all 
r > rg by the ratio test. The proper distance now be- 
comes 



£{r) 



kro 



,b['nr-^T ^midr. 



(19) 



To study the behavior of £{r) as a function of p, we let 
rg = 4, /c = 2, a = 1, and b = -0.20. To show how 
£{r) drifts away from the original values (corresponding 
to p = 1), l/£{r) is plotted against p, as shown in Fig. 1. 
The graph resembles the right half of a normal curve. So 



B. Other solutions 

One of the main conclusions in Ref. ^14] is the extreme 
difficulty in obtaining exactly solvable wormhole models 
without getting an unwanted event horizon. While the 
question of exact solutions does not arise in this study, 
based on Eq. and numerical trials, choosing an "arbi- 
trary" $'(r) would not ordinarily yield an event horizon. 
The choices for $ are therefore enormous. 

Another fortunate circumstance is the occurrance of 
A(r) = a + b\n{r — rg) as a best-fitting curve. This 
function causes e^^''-' to go to infinity relatively slowly, 
producing a small proper distance, as exemplified by 
Eq. So for small |$'(?')| we are dealing with a 

situation that can hardly be improved. As |<l''(r)| in- 
creases, however, the In-functions used so far produce an 
ever poorer fit. The question is whether this deviation 
from the In-form can be quantified. A related question 
might be: given the abundance of solutions and the re- 
sulting possibility of naturally occurring wormholes, can 
one estimate the likelihood of finding one that is actually 
traversable by humanoid travelers? 




p-values 



FIG. 1: A measure of traversability. 

if we regard l/i{r) as a measure of traversability, then, 
according to Fig. 1, the likelihood of encountering a nat- 
urally occurring GCG wormhole that is also traversable 
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may be higher than one should expect, perhaps even as 
high as 25% (corresponding roughly to 1 > p > 0.6). 

This outcome is fairly typical: increasing the throat 
size ro has little effect on the shape of the graph, while 
the values of k and a have no effect at all. As b decreases 
toward —0.5, the graph falls off more rapidly, but the per- 
centage of traversable wormholes appears to remain be- 
tween 10% and 20%. Fortunately, even the worst-fitting 
cases have led to values between —0.35 and 0. 

So far we have assumed that a in the equation of state 
p = —Ajp'^ is close to unity, thereby producing the best- 
fitting curves. As a decreases, one sees the same kind of 
falling off in the measure of traversability that appears 
in Fig. 1. So the best chance for obtaining a traversable 
wormhole is a GCG that is close to the original Chap- 
lygin gas. However, this falling-off behavior can be par- 
tially compensated for by decreasing the constant A in 
the equation of state. This can also be seen from Eq. 
an increasing 1 j a is compensated for by decreasing A. So 
a traversable wormhole could in principle exist for small 
a, provided that the size of A can be controlled (for exam- 
ple, by an advanced civilization). In that case, however, 
the inequality ([5]), p < becomes a constraint 

on the wormhole material. 



C. Junction to an external spacetime 

Observe that any curve of the form A(r) — a-|-51n(r — rg) 
with a > and 6 < will eventually cross the horizontal 
axis. So our spacetime is not asymptotically flat, the 
same problem that occurs in Refs. [13 and [ii]. We 
would like to join A(r) smoothly to a positive curve that 
goes to zero, as in Ref. [i8|. 

Assume that the extended curve has the form g{r) — 
K/r"^, starting at some r = ri. We require that g{ri) — 
A(ri) and ^'(ri) = A'(ri). Eliminating K and n, we 
obtain 



g(r) = A(ri) 




riA'(ri)/A(ri) 



Since we are concerned mainly with the region around 
the throat, let us opt for an early cut-off, say at r = 2.2. 
(Besides, according to Ref. [l^, r lies in a fairly nar- 
row range, which is also apparent from the numeri- 
cal/graphical output.) As an example, for the case 
$'(r) = 0, we have A(r) = 1.20 - 0.23 ln(r - 2). Then for 
r = 2.2, 

/ r \ -1-61 
9ir) ^ 1.57 (-) . 

Joining A(r) to a positive curve going to zero makes 
the spacetime asymptotically flat. 



V. COMPARISON TO LOBO'S SOLUTION 

The case $'(r;) = allows a comparison to Lobo's ex- 
act solution tia]. From A(r) = 1.20 - 0.23 ln(r - 2), we 
calculate 

6i(r)=r[l-e-2-4(r-- 2)0-46]. 
Lobo obtains 

rfi4 Att'^ 1 
. -J ^0 

On the above interval [2,2.2], using the same A and rg, 
bi{r) = b2{r) after rounding off to two significant figures. 

VI. SUMMARY 

This paper extends the general analysis in Ref. 
to wormholes supported by generalized Chaplygin gas 
(GCG). The function A(r) in line element ([2|) is obtained 
numerically from the Einstein field equations and the 
equation of state. The conclusions are: 

1. The wormhole spacetime meets the flare-out condi- 
tions at the throat. 

2. Qualitatively, A(r) remains the same for a wide 
range of redshift functions <i>(r). Compared to the 
phantom energy supported wormholes in Ref. p^ . 
event horizons are much less likely to occur. 

3. Redshift functions with relatively small [^''(r)! 
yield regression equations closely fitting the form 
A(r) = a + 61n(r — tq). 

4. Item 3 implies the existence of macroscopic worm- 
holes with low tidal forces and relatively short 
proper distances near the throat, making the worm- 
holes traversable by humanoid travelers. (This cri- 
terion does not include proper traversal times.) 

5. The junction to an external solution produces a 
spacetime that is asymptotically flat. 

6. It is assumed above that a is close to unity in the 
equation of state p — —A/p°'. A smaller a produces 
less favorable results but can be compensated for by 
a smaller A. However, a smaller A would tighten 
the constraint p < j4^/("+^' on the wormhole. 

7. The abundance of solutions suggests that GCG 
wormholes may occur naturally. A rough measure 
of traversability implies that the chances of finding 
one that is traversable by humanoid travelers may 
be fairly good. 
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